Photon detection by current-carrying superconducting film: A time-dependent 

Ginzburg-Landau approach 
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We study dynamics of the order parameter in a superconducting film with transport current after 
absorption of a single photon. The system from the time-dependent Ginszburg-Landau equation, 
Poisson's equation for an electrical potential and the heat diffusion equation was solved numerically. 
For each photon energy in the absence of fluctuations, there exists a corresponding threshold current 
below which the superconducting state is stable and no voltage appears between the ends of the 
film. At larger currents the superconducting state collapses starting from the appearance of a 
vortex-antivortex pair in the center of the region with suppressed superconducting order parameter 
which has been created by the absorbed photon. Lorentz force causes motion of these vortices that 
heats the film locally and gives rise to a normal domain. When biased with the fixed current, the 
film latches in the normal state. In the regime when current via superconductor may change, which 
is more relevant for experiments, the normal domain exists only for a short time resulting in the 
voltage pulse with the duration controlled by the kinetic inductance of the superconducting film. 

PACS numbers: 74.25.0p, 74.20.De, 73.23.-b 



I. INTRODUCTION 

Despite the large number of both experimental (see for 
example iMZl and references therein) and theoretical 0- 
[Tl| works on the superconductive single-photon detectors 
(SSPD) still there are some questions about mechanism 
of photon detection by superconductive film carrying the 
transport current. The original understanding of the de- 
tection mechanism is following - after absorption of the 
single photon the hot spot is formed in the superconduct- 
ing film. It locally destroys superconductivity and leads 
to concentration of the current density outside the hot 
spot due to decrease of the effective width of the film Q . 
If transport current is close to the depairing current then 
the current density outside the hot spot area may exceed 
the depairing current density and the superconducting 
state becomes unstable, leading to the voltage response. 

The quantitative analysis of initial stage of the hot spot 
formation in the existing theoretical models is based on 
the solution of the diffusion equation for nonequilibrium 
quasiparticles jst-TlOl]. It was postulated that when the 
number of the nonequilibrium quasiparticles in the hot 
spot exceeds some critical value this region could be con- 
sidered as a normal one and the superconducting current 
is forced to flow around it. In the refined model @ it 
was supposed that even partial suppression of the su- 
perconducting order parameter in the hot spot leads to 
current enhancement outside that region and to insta- 
bility of the superconducting state and formation of the 
normal domain. The further evolution of the normal do- 
main is usually studied by using the heat diffusion equa- 
tion for effective temperature of quasiparticles coupled 
with equation describing the embedding circuit |lll4l3l |. 
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In the majority of previous models it was implicitly as- 
sumed that the magnitude of the superconducting order 
parameter | A| changes instantly in time if local tempera- 
ture T(r, t) > Tc or local superconducting current density 
j{r,t) exceeds depairing current density jdep- But it is 
well know that |A| has finite relaxation time r|A| and 
in some cases T|a| could be comparable with electron- 
phonon inelastic relaxation time Te-ph • Because en- 
ergy relaxation of nonequilibrium quasiparticles in the 
hot spot occurs on the same time scale (or even much 
shorter due to diffusion of the quasiparticles) it is clear 
that finite r|A| 7^ should influence the photon detection 
process [Tsj . 

Another interesting and unresolved question is what 
kind of instability of the superconducting state occurs 
due to appearance of the hot spot region. Is it gradual 
suppression of the order parameter outside the hot spot 
due to current concentration [l^ or nucleation of the 
vortex-antivortex pair inside the hot spot (T6|? 

In our work we use the simplest approach where ef- 
fect of finite T|a| is taken into account and stability of 
the superconducting state is analyzed self-consistently. 
The dynamics of the superconducting order parameter 
is studied on the basis of the time-dependent Ginzburg- 
Landau equation. This equation is coupled with the heat 
diffusion equation for the effective temperature of the 
quasiparticles and Poisson's equation for the electrical 
potential. We consider the current bias regime as well as 
the regime when current via superconductor may change 
due to presence of the shunt resistance and take into ac- 
count the finite kinetic inductance of the film. 

Within this model we show that incoming photon cre- 
ates the finite size region with partially suppressed or- 
der parameter. We find that even for infinite supercon- 
ducting film such a state becomes unstable without any 
fluctuations with respect to appearance of the vortex- 
antivortex pair (or single vortex if photon is absorbed on 
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the edge of the film) at threshold current less than depair- 
ing current. Motion of the vortex and antivortex in oppo- 
site directions under the Lorentz force heats the sample 
(if the threshold current is not too small). As a result, 
the normal domain appears which either expands over 
the whole film (current bias regime) or shrinks and dis- 
appears (when current via superconductor may change) 
resulting in the voltage pulse. Our result supports the hy- 
pothesis of Ref. [l^ that the single photon can nucleate 
the vortex-antivortex pair in current-carrying supercon- 
ductor and their motion provides the voltage pulse which 
could be detected. Our model also confirms the exper- 
imentally observed smeared red boundary in the single- 
photon detection. 

The paper is organized as follows. In section II we 
present the theoretical model. The results of the nu- 
merical calculations and simple analytical estimations are 
given in Section III. In section IV we discuss the relation 
of our results with an experiment and in section V we 
present our conclusions. 



II. MODEL 

In our work we do not study the initial part of the 
detection process when the single photon is absorbed by 
the electron. We use approach of the effective temper- 
ature [l7ji, which is valid when the thermalization time 
(which is proportional to the electron-electron inelastic 
relaxation time Te_e) is shorter than the inelastic relax- 
ation time due to electron-phonon interactions T^-ph ■ We 
assume that during initial time interval ~ Te_e after ab- 
sorption of the photon the electron-electron interactions 
creates hot spot with radius Rinit ^ -^e-e — (^'^e-e)^^^ 
(D is a diffusion constant) and with local temperature 
To + AT (To is a bath temperature) where AT is deter- 
mined from the energy conservation 



27rhc/X = ATirRj^.^dCy 



(1) 



Here A is a wavelength of the electromagnetic radia- 
tion, h is a Planck constant, c is a speed of light, d is 
a thickness of the film and is a heat capacity of the 
quasiparticles (for simplicity we take Cy as in the normal 
state at T = Tc). 

Time and space evolution of the temperature in the su- 
perconducting film we find from the heat diffusion equa- 
tion 



dT 
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Cy 
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where p„ is a normal state resistivity, j„ — —'^^jpn is 
a normal current density and is a electrostatic poten- 
tial. Here we assume that the phonons are in equilibrium 
with the bath and energy relaxation occurs due to inter- 
action with phonons. Our calculations show, that ini- 
tial destruction of superconductivity occurs on timescale 
shorter than Te-ph and therefore at initial stage of dy- 
namical response of |A| one may neglect the heating 



of phonons (in our model we neglect possibility of the 
phonon heating during initial t < r^-e stage of hot spot 
formation) . 

To study the dynamics of the order parameter A = 
lAje**^ we use the time-dependent Ginzburg-Landau 
equation 
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dt 



A = eGL(0)- 
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A, (3) 



where ^gl(O) = [nhD/SkBT^^/^ and Agl(O) = 
4:kBTcU^/'^/Tr {u ~ 5.79 - see ^) are the zero tem- 
perature Ginzburg-Landau coherence length and the or- 
der parameter respectively. Characteristic time relax- 
ation of the order parameter described by Eq. (3) is 
T|A| = nh/8kB{Tc - T). Although Eq. (3) is quanti- 
tatively valid only near critical temperature of the su- 
perconductor (at T > 0.9Tc when re_e <C Te-ph and 
Te_e ^ '''|A|) we use it to model the dynamics of the su- 
perconducting condensate at lower temperatures to find 
some qualitative results. 

We should complete Eqs. (2-3) by equation for the 
electric potential which comes from the conservation 
of the full current div{js + jn) = 



A(p = pndiv{js) 
here j, = Imag(A*VA)/(4e/cBT,p„). 



(4) 
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FIG. 1: The model geometry: the superconductive film is 
placed between two normal bulk contacts; (a) - the photon is 
absorbed in the center of the film, (b) - the photon is absorbed 
on the edge of the film. 

To model the response of the superconducting film af- 
ter absorption of the single photon we consider the model 
geometry which is present in Fig. 1. We need the normal 
contacts (A = 0, d^p/dx — —pnl/{wd) a.t x — ±1/2) to 
inject the current to the superconducting film in our nu- 
merical calculations and which are kept at the bath tem- 
perature To (T|^^-|-;/2 = Tq). The current and heat do 
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not flow through the lateral edges of the film {dT/dy = 0, 
dip/dy = 0, dA/dy = at y = ±u;/2). To neglect the 
influence of the N-S boundaries (for example the mo- 
tion of the NS boundary) on the dynamical processes in 
the superconducting film we artificially enhance the su- 
perconducting order parameter in the regions marked by 
dark blue color in Fig. 1 by introducing locally higher Tc 
(the width of these regions is larger than the penetration 
depth of the electric field from the normal contact and is 
equal to 5Cgl(0)). 
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FIG. 2: The equivalent scheme of the superconducting detec- 
tor. The superconductor is modelled by kinetic inductance 
Lk and resistance Ra which appeared due to absorbing the 
photon. The shunt has resistance Rshunt- 

To model real experiments we consider electrical 
scheme which is shown in Fig. 2. Here Lj- is the kinetic 
inductance of the superconducting film, Rs corresponds 
to the resistance of the superconductor in the resistive 
state (in the model geometry - see Fig. 1) and Rshunt is 
the shunting resistance. For this case we have to find cur- 
rent Is which flows via superconductor from the solution 
of the following equation 



L^dls_ _ . 
dt 



Vs 



(5) 



where the voltage drop over superconductor Vs (over the 
blue region in Fig. 1) should be found from the solution 
of Eq. (4) with boundary condition dip/dx\x=±L/2 = 

~PnIsl{wd). 

In numerical calculations we use the dimensionless 
units. The order parameter is scaled in units of AGi(O), 
temperature is in units of Tc and coordinate is in units 
of ^gl(O). Time is scaled in units of ro = irh/SkBTcU, 
electrostatic potential is in units of tpo = h/2eTQ and cur- 
rent density is in units of jo = V2ep„ro<^GL(0) (depairing 
current density in these units is jdep/jo — (4/27)^/^(1 — 
T/T,)3/2). 

To solve Eqs. (2,3,5) numerically we use Euler method 
and to solve Eq. (4) - Fourier analysis and cyclic reduc- 
tion method. In numerical calculations we flrst apply the 
finite current and wait until all relaxation processes con- 
nected with the current induced suppression of the order 



FIG. 3: Time dependence of the magnitude of the order pa- 
rameter (a) and temperature (b) in the center of the hot spot 
(which coincide with a center of the film) for two values of the 
transport current / = 0.89/dep and / = 0.9/dep. The width of 
the film w — 525gl(0), the local initial increase of the tem- 
perature AT — 2.3Tc (A ~ 6.5/im). In the inset we show 
contour plots of the magnitude of the order parameter in the 
film at different times marked by the numbers on the black 
solid curve. 



parameter stops. Than at some moment of time we in- 
stantly increase the temperature by AT in the circle or 
semicircle area inside the superconductor (see Fig. 1) and 
studied the dynamical response of the system. The pa- 
rameters of the film are: length / = 60^gl(0) and width 
w is varied from 13^gl{0) up to 78 ^gl(O). 



In our calculations we use parameters typical for NbN 
SSPD [|: a = 2AmJcm-^K-^ , r^-e = 7ps, D = 
0.45cmVs, Cgl(O) = 7.5nm, = lOK, r^-ph = 17 ps. 
At these parameters Le_e ~ ISnm and tq = 0.052ps. In 
test calculations we consider two values for Rinit = 18nm 
and Rinit = 9nm which are close to Le-e and found that 
the results (in particular the value of the threshold when 
the voltage appears) differ only slightly. The presented 
below results are obtained with Rinit — 9nm. For this 
radius and thickness of the film d = Anm the range of 
AT = 0.3 — 12.8Tc corresponds to the wavelengths A ~ 
1.3 — 50/im. The bath temperature To is equal to Tc/2. 
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III. RESULTS 
A. Regime with constant current 

At first we consider current bias regime (when Is = I 
and Ishunt = in Fig. 2 for Rshunt oo). In Fig. 3 
we present time dependence of the magnitude of the or- 
der parameter and effective temperature of quasiparticles 
in the center of the film with width w = 52^ci(0) and 
for the situation depicted in Fig. 1(a) for two close val- 
ues of the transport current. Notice that suppression of 
the order parameter in the center of the hot spot needs 
finite time (see Fig. 3(a)). During this time the local 
temperature in the center of hot spot decreases (see Fig. 
3(b)) due to diffusion of the nonequilibrium quasiparti- 
cles and energy transfer to phonons. When the current 
is smaller than the threshold value (we call it the de- 
tecting current Id) the order parameter after reaching 
some minimal value starts to grow. In this case the time 
averaged voltage response is zero. The larger current de- 
stroys superconducting state. In this case |A| oscillates 
in the center of the hot spot with the amplitude which 
decays in time. Each oscillation of |A| corresponds to 
nucleation of one vortex-antivortex pair. Motion of the 
vortex/antivortex in opposite directions (see inset in Fig. 
3(a)) heats the superconductor via Joule dissipation and 
the local temperature increases. It results in appearance 
of the growing resistive domain (see inset in Fig. 3(a)) in 
the regime of the constant current at chosen parameters. 

In Figs. 4-5 we show the time evolution of the order 
parameter in the films with smaller width 'At I > Id and 
one can see qualitatively the same scenario of the order 
parameter dynamics. The dependence of the detecting 
current on AT (i.e. on the energy of the absorbed pho- 
ton) and width of the film is shown in Fig. 6. The 
detecting current decreases with the increase of the pho- 
ton energy and its value depends on the position where 
the photon is absorbed (in the center or on the edge of 
the film). For fixed photon energy the ratio Id/Idep ini- 
tially grows with increasing width of the film and than 
saturates for large w (see Fig. 6(b)). 

We shall note that for high energy photons (large 
AT) and relatively narrow film the detecting current is 
much smaller than depairing current (see Fig. 6(a)). In 
this case the Joule dissipation could be weak, the nor- 
mal domain does not appear and superconductivity re- 
covers after nucleation of several vortex-antivortex pairs 
in the hot spot area. For film with w = 13^g'l(0) 
and AT = 12.8Tc the normal domain appears only at 
/ > 0.36Idep {Id — 0.23/dep - see Fig. 6(a)) which is close 
to the value of the current when the heat dissipation and 
heat removal are equal to each other 

PnJheat _ — Tp 

and for our choice of parameters Iheat = j heated ~ 

0.23/dep- 




FIG. 4: Time dependence of the magnitude of the order pa- 
rameter in the center of the hot spot (which coincides with 
the center of the film). The width of the film w = 26£,gl{0), 
the bias current 1=0.86 Idep, the local initial increase of the 
temperature AT — 2.3Tc (A ~ 6.5/im). The inset show con- 
tour plots of the magnitude of the order parameter in the film 
at different times marked by the numbers on the black solid 
curve. 
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FIG. 5: Time dependence of the magnitude of the order pa- 
rameter in the center of the hot spot (which coincides with 
the center of the film). The width of the film w = 13^gl(0), 
the bias current I = 0.75/dep, the local initial increase of the 
temperature AT — 2.3Tc (A ~ 6.5/im). The inset show con- 
tour plots of the magnitude of the order parameter in the film 
at different times marked by the numbers on the black solid 
curve. 



According to our numerical calculations the voltage 
response appears when the vortex-antivortex pair is nu- 
cleated in the center of the hot spot. To get insight 
why it occurs we consider the following simple model. 
Let us model the region with suppressed |A| in the hot 
spot by circle of radius R and we assume that |A| is 
spatially uniform and has value Ai„ inside the circle 
and in the rest of the infinite thin superconducting film 
|A| = AoMt > Ai„. We are interested how transport cur- 
rent is distributed in such a superconducting system and 
when superconducting vortex free state becomes unsta- 
ble. For simplicity we neglect the proximity effect (which 
is reasonable when i? 3> ^) and use the London model 
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FIG. 6: (a) The dependence of the detecting current on the 
instant increase of the temperature in the circle with radius 
Rinit for narrow and wide films with two positions of the 
photon absorbtion (in the center and on the edge of the film), 
(b) The dependence of the detecting current on the width of 
the superconductive film for three values of AT/Tc=0.8, 2.3, 
3.8, corresponding to three wavelengths of the electromag- 
netic radiation A =18.8, 6.5 and 3.9 [im respectively (photon 
is absorbed in the center of the film). 



is — I ApV0/(4eA:BTcp,i). Distribution of the supercon- 
ducting current can be found from the current conser- 
vation divjs = 0. As a result we obtain that inside the 
circle the supervelocity v ~ V0 is larger than v in infinity 



2w„ 



1 + 7' 



(7) 



(7 = Ai„/Aout) and it is locally enhanced outside the 
circle 



1 



d2 1 2 
/t 1 — 7 



1 



r 



,r>R 



(8) 



where the distance is measured from the center of the 
circle and we present result at angle a = tt/2 between 
the direction of the current and radial vector in polar 
system of coordinate. 

One may find corrections to Eqs. (7,8) for film with 
finite width in the limit when 2R/w <C 1 and the circle 
is placed in the center of the film (see Fig. 1(a)). As- 
sume that for finite film with 2R/W < 1 the Eqs. (7,8) 



are approximately valid, but coefficient Voo we replace by 
unknown v* which we find from the conservation of the 
full current 



/ Al^^Voowd = 2d 



1 + 7' 



-dy 



(9) 



,w/2 
-2d / a2 

Jr 



1 + T 



As a result we find 

V* = Voo/ 1 




(10) 



Note that Eq. (10) is also valid (with replacement 
2R/w — R/w) for the case when semicircle of radius 
R with suppressed | A| = Ai„ is placed on the edge of the 
film (see Fig. 1(b)). We have to stress that coefficient 
in front of term {2R/w)'^ in Eq. (10) is approximately 
valid (up to coefficient of order of unity) and correct value 
should be found from the expansion of the exact result 
in series with small parameter 2R/w. 

Because Ai„ < Ao„t and Vm > Vout the superconduct- 
ing Meissner (vortex free) state first becomes unstable 
inside the circle(semicircle). Using value of the critical 
supervelocity Vc ^ |A| for instability of spatially uniform 
superconducting state which follows from stationary Eq. 
(3) we find (by equalizing Vm = Vc and using Eq. (10)) 



Idep 



7(1+7') 



1 - 



2R 

w 



1-7^ 
1+7' 



(11) 



The current Ipair is our estimation for the current when 
the vortex-antivortex pair is nucleated inside the circle 
(with replacement 2R/'w — > R/w it corresponds to the 
threshold current when single vortex is nucleated in the 
center of the semicircle on the edge of the film). It is 
worth to mention here that it is not the current when the 
resistive state appears in the sample, because to escape 
the circle(semicircle) the vortex and antivortex should 
overcome the energy barrier connected with jump in |A|. 
To estimate this critical current we assume that the vor- 
tices may leave the circle(semicircle) when the superve- 
locity (averaged over finite region ^ £,{T) near the edge 
of the circle) is equal to Vc- Using Eqs. (8, 10) it is easy 
to find that 



^res 
Idep 



= 1- 



2R 



1 



7 



1 + 7' 



/ 1 + 



R 



1-7' 



One can see that /„, 



R + C(T) 1 + 72 

(12) 

< Ires (they are equal when 
Am = Ao„t and 7=1). Both critical currents decrease 
with decreasing Ai„ and Ipair = and Ires = Idep{^ — 
AR^/w'^)/2 when R > ^(T) and 7 = (it corresponds to 
the normal state of the circle) . 

To complete the analytical analysis we have to corre- 
late the radius of the region with suppressed | A| and Aj„ 
with energy of the incoming photon. Let us assume that 
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the spatial and time dependence of the temperature after 
photon absorbtion is described by the following expres- 



sion 



T{r,t) 



AnDt 



Tn 



(13) 



which is solution of the Eq. (2) with replacement of heat- 
ing term by f35{t)5{y') {(i — 2Trhc/ XCyd) which describes 
the energy delivered by photon to the quasiparticles at 
the moment t=0 and in the point r=0 (we also neglect 
the last term in Eq. (2) because we are interested in time 
interval of about T|a| <C T^-ph after photon absorption). 

Local enhancement of temperature leads to suppres- 
sion of the order parameter in the hot spot. We may 
estimate it by using Eq. (3) where we neglect for sim- 
plicity the term with the second derivative 



^|A|(0) 



d\A[ 
dt 



T 



lAI 



(0)^ 



(14) 



One can see that while left hand side of Eq. (14) is 
negative the order parameter decreases. Because T is 

maximal in the center of hot spot and decreases in time 
it is reasonable to suppose that the order parameter stops 
to decrease when T = Tc in the center of the hot spot. 
Using Eq. (13) we find that it occurs at 



5t = 



47rD(Te - To) T, 



(15) 



where we used Eq. (1) to express /3 = ATnRf^^^ 
and Rinit = 1-2^gl(0) via parameters of our numerical 
model. Using this result and Eq. (13) we may estimate 
the size of the region where the order parameter is sup- 
pressed 



/3 



7r(Tc - To] 



AT 



(16) 



increase of the photon energy R increases and 7 decreases 
providing decrease of L^es and Ipair (compare with Fig. 
6(a)). 

Secondly, above analytical results also explain decrease 
of detecting current with decreasing width of the sample 
(compare Fig. 6 and Eq. (12)). Thirdly, when the pho- 
ton is absorbed on the edge of the film it creates the semi- 
circle with larger radius R' = \/2R (in comparison with 
photon absorbed in the center). It results in larger de- 
tecting current than for photon absorbed in the center 
of the film (see Eqs. (11,12) with replacement 2R/w by 
R' /w and for high energy photons when 7 <Si 1). Note 
that the effect is stronger for films with smaller width 
(compare with Fig. 6(a)). 

For photons of relatively small energy {AT/Tc > 1) 
which create the circle (semicircle) with small effective ra- 
dius R ~ ^(T) and wide film {w ^ ^{T)) the situation is 
different. In this limit the correction factor due to finite 
w in Eqs. (11,12) is small and one can see that Id is 
smaller for the photon absorbed on the edge than for the 
photon absorbed in the center of the film (due to differ- 
ence in radiuses and in 7 which is finite). It correlates 
with our numerical results for low energy photons (small 
AT) and wide film (see Fig. 6(a)). 

We have to note that because of temperature gradient 
and proximity effect the distribution of the order param- 
eter is nonuniform in the hot spot formed by the pho- 
ton. This factor was not taken into account in above 
model and it brings the quantitative difference between 
our numerical and analytical results. Our numerical cal- 
culations show that at any considered photon energy the 
order parameter is finite in the hot spot at the moment 
when the first vortex-antivortex pair is nucleated. After 
nucleation the vortex and antivortex becomes immedi- 
ately unbound and can move freely across the supercon- 
ducting film. Therefore found above values for Ipair and 
Ires could be considered as low and upper thresholds for 
true detecting current. 



From Eq. (13) one may find suppression of |A| in hot 
spot at r < ii by moment t = St. Below we assume that 
Air, = \A\{r = ^,t = St). Using Eqs. (13,14) we find 
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Aout 

~ exp 



exp 



i_ r'' 

r\A\{0)Jo 



-In 



dt (17) 



4<GL(0)2rc \ttS,Gl{Q)''TcJ J 

which is approximately valid for photons with 

(/3MGL(o)2r,c. Ar/r,>i). 

For photon absorbed on the edge (which creates the 
semicircle) above results are also valid with the replace- 
ment of P hy 2(3. Finite width of the film does not affect 
above results too much while w ^ 2R (or w ^ R for the 
semicircle) due to exponential decay of temperature at 
r>R. 

Combination of Eqs. (16,17) and Eqs. (11,12) quali- 
tatively explains our numerical results. First of all with 



B. Regime with changing current 

In the experiments the circuity shown in Fig. 2 is usu- 
ally used to prevent the latching of the superconductor 
in the normal state. As a result with appearance of the 
voltage drop over superconductor the current via the su- 
perconducting sample decreases and it switches back to 
the superconducting state. 

In our calculations we use Rshunt=5Q Ohm. The ki- 
netic inductance is evaluated using the following expres- 
sion: 



Lik = 3— 

wa 



(18) 



where Al is a London penetration depth. 

In Fig. 7(a) we present the time dependence of the 

voltage drop Vg via the superconductor (resistance Rg 
in Fig. 2) for two films with lengths 1=100 /xm and 
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FIG. 7: (a) The time dependence of the voltage drop via 
the superconductive film calculated for two lengths: 1=100 
/xm and 1=200 ^m. The inset shows the distribution of the 
magnitude of the order parameter along the film (at y=0) at 
different moments in time, (b) The time dependence of the 
temperature in the center of the film (hot spot) for the same 
films. The inset shows the distribution of the temperature 
along the film (at y=0) at different moments in time. The 
width of the film w = 13^gl(0), the bias current / = O.G/dep, 
the local initial increase of the temperature AT = S.STc (A ~ 
3.9/^m). 

200 /iTO (w=100 nm, d=4 nm and Al = 400 nm 
It is seen that with decrease of the kinetic inductance 
the duration and the amplitude of the voltage pulse be- 
comes shorter and smaller correspondingly. The reason 
is simple - for smaller Lk current via superconductor de- 
creases faster, the temperature inside the normal domain 
increases slower (see Fig. 7(b)) and it takes less time 
to cool the sample up to bath temperature Tq (see Fig. 
7(b)). 

For the longest film the duration of the voltage pulse 
is about 900 tq which is much larger than the typi- 
cal time interval between consequent nucleation of the 
vortex- antivortex pair At ~ IOtq (we roughly estimate it 
as a time interval between nucleation of the second and 
third vortex-antivortex pairs - see Fig. 5). Therefore at 
least 90 vortex-antivortex pairs are nucleated during this 
voltage pulse. 

In Fig. 8(a) we present the time dependence of Vs for 
photons of different energy and position where it is ab- 
sorbed. Notice that the shape of the voltage pulse and 
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FIG. 8: (a) The time dependence of the voltage via the su- 
perconductive film with length 1=200 and width w = 
26^Gi(0) for photons of different energy (different AT) and 
position of the absorbtion (in the center or on the edge of 
the film), (b) The time dependence of the current via shunt 
resistance (the bias current / = 0.887dep)- 



current Ishunt (see Fig. 8(b)) via shunt resistance (which 
is measured in the experiments with SSPD) slightly de- 
pend both on AT and absorption position. At times 
larger than IOOOtq when the voltage drop via supercon- 
ductor is equal to zero (see Fig. 8(a)) the current Ishunt 
decays with characteristic time r = Lk/ Rshuntc^ — 
5 • 10'*To for our choice of parameters. 

We have to stress that our results for voltage pulse 
via superconductor (duration of the pulse, its ampli- 
tude) are obtained with neglecting the back action of 
the 'heated' phonons on the electrons. For our parame- 
ters {Te-ph — 327ro) the duration of the pulse is about 
"^Te-ph- But heating of the phonon subsystem occurs 
over time t 3> T^-ph '203. Therefore we expect that hot 
phonons do not strongly influence the duration of the 
voltage pulse and its amplitude. And it certainly should 
not influence the value of the detecting current because 
it is determined from the stability of the superconduct- 
ing state with respect to appearance of the single vortex 
or vortex-antivortex pair during evolution of the photon 
induced hot spot. 
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IV. RELATION TO AN EXPERIMENT 

The model which we use (Eqs. (1-4)) is strongly over- 
simplified. First of all it does not take into account lose 
of the energy of photon at initial stage of formation of 
the hot spot on time scale r^-e- Secondly we neglect 
temperature dependence of C^, oversimplify the energy 
transfer to the phonons (last term in Eq. (2)) and we 
did not take into account possibility of direct partial de- 
struction of the superconducting order parameter by in- 
coming photon (in our approach |A| is influenced only 
via effective temperature of quasiparticles). Moreover, 
the time-dependent Ginzburg-Landau equation (Eq. (3)) 
is not quantitatively correct at temperatures lower than 

O.OTc and (within quasiequilibrium approach) when 
Te-e ^ ''"|A|- Therefore direct quantitative comparison 
of our results with experiment (at least at low tempera- 
tures) looks speculative. 

But despite this we believe that the used model catches 
the main physical mechanism of photon detection by 
current-carrying superconducting film which is the fol- 
lowing. The incoming photon partially suppresses order 
parameter in the finite region. It leads to redistribution 
of the current density (supervelocity) in the film and such 
a state becomes unstable [without any fluctuations and 
if current is large enough but smaller than depairing cur- 
rent) with respect of appearance of the unbound vortex- 
antivortex pair (if photon is absorbed far from edges) or 
single vortex (if photon is absorbed near the edge of the 
film). Lorentz force causes motion of these vortices that 
heats the film locally and gives rise to a voltage pulse. 

We hope that our results could be used for understand- 
ing (not for direct fitting) of some experimental results. 
For example found dependence of the detection current 
on the position of the hot spot (see Fig. 6(a) and our an- 
alytical model) may be used for qualitative explanation 
of the monotonous decrease of the detection efficiency 
(DE) with the decrease of the energy of the incoming 
photon 0,111, HI]- Indeed, if we fix the current (for ex- 
ample on the level / — Idepf^ - see Fig. 6(a)) and start 
to decrease the energy of the photon, the edge region of 
the superconducting film first is "switched off' from the 
detection process (when AT ~ 8Tc which corresponds to 
A ~ 1.9fim - see Fig. 6(a)) and at AT ~ ST^ (A ~ 2.9^m) 
the central region of the superconducting film stops to 
detect photons. As a result there is finite range of the 
wavelengths AA ~ l/um at (/ = Jdep/2) where the de- 
tection efficiency gradually changes (qualitatively such a 
behavior was observed in Refs. 0, [Ij, [13 )■ 

If we fix the energy of photon then with increase of 
the current the central region of the film first starts to 
detect the photons and than the edge regions join the de- 
tection process. Therefore there is a finite interval of the 
currents {51) within which DE gradually grows with cur- 
rent increase (qualitatively such a behavior was observed 
in Ref. [25|). The real samples have different kinds of 
imperfections (variations of the thickness, width, mate- 
rial parameters, turns) having their own values of the 



detecting current and it obviously affects 51 at low cur- 
rents. Our result show that even in ideal samples with 
no imperfections and at low temperatures (when effect of 
fiuctuations is rather small) there will be finite 51 con- 
nected with presence of the edges. And our prediction is 
that the wider the sample the narrower is this interval of 
currents (see Fig. 6(a) and Eq. (12)). 



V. CONCLUSION 

In our work wc use the quasiequilibrium approach and 
describe the deviation from the equilibrium in terms of 
the effective temperature of the quasiparticles which de- 
pends on time and coordinate. We assume that the ab- 
sorbed photon creates initially the hot spot in the su- 
perconducting film with radius Rinit and local enhance- 
ment of the quasipartcile temperature by AT which is 
proportional to the energy of the photon. The tempo- 
ral and spatial evolution of the effective temperature and 
superconducting order parameter in the superconductor 
we study by numerical solution of the time-dependent 
Ginzburg-Landau equation coupled with Poisson's equa- 
tion for an electrical potential and heat diffusion equa- 
tion. 

We show that for photon of fixed energy the voltage 
response appears only at the current larger some critical 
one (we call it as detecting current Id) when the vortex- 
antivortex pair is nucleated in the center of the photon 
induced hot spot. Motion of the vortex and antivortex 
in opposite directions heats the superconductor and leads 
to the appearance of the growing normal domain when 
detecting current Id is about of depairing current. For 
high energy photon and narrow film Id Idep the mo- 
tion of the vortices-antivortices does not heat the super- 
conductor and sample goes back to the superconducting 
state after nucleation of several vortex-antivortex pairs 
at I Id even in regime of constant current. 

We find numerically that with increasing the width of 
the superconducting film Id increases and stays less than 
Idep even for infinite film. We also find that the detecting 
current for photon absorbed on the edge of the film differs 
from Id for photon absorbed in the center of the film. 

We develop simple analytical model to explain above 
results. We assume that absorbed photon creates in the 
superconducting film of finite width the region (in the 
form of circle or semicircle) where the superconducting 
order parameter |A| is partially suppressed and in frame- 
work of the London model we study the current redistri- 
bution and stability of the current-carrying state in such 
a system. We find that the superconducting vortex free 
state becomes unstable in the region with suppressed su- 
perconductivity at Ipair < Idep and resistive state ap- 
pears at larger current Ipair < Ires < Idep- Our ana- 
lytical model predicts different values for the detecting 
current of the photon absorbed on the edge and in the 
center of the film in agreement with our numerical results. 
The last effect is connected with different redistribution 
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of the supervelocity (current density) inside and outside 
the circle placed in the center of the film and the semicir- 
cle placed on the edge of the same film. We expect that 
our analytical results within the London model are valid 
for arbitrary temperatures contrary to the results based 
on the time-dependent Ginzburg-Landau and heat con- 
ductance equations which are strictly valid at T > 0.9Tc 
and when Tg-e < T|a| and Tg-e < r^-ph- 

To model operation of real superconducting single pho- 
ton detector we consider the scheme with the resistance 
which is switched on in parallel to the detector and take 
into account the large kinetic inductance of real SSPD. 
We find that duration and amplitude of voltage pulse de- 
crease with the decrease of the kinetic inductance while 
the detecting current practically does not change. We 



also find that the shape of the voltage pulse weakly de- 
pends on the energy of the absorbed photon and on the 
place of the absorbtion for the homogeneous film. 
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